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ABSTRACT 



> 

Q\ ' We present the results of a second epoch of 18 cm global Very Long-Baseline Inter- 

£NJ , ferometry (VLBI) observations, taken on 23 February 2001, of the central kiloparsec 

^sO ■ of the nearby star burst galaxy Messier 82. These observations further investigate the 

CO structural and flux evolution of the most compact radio sources in the central region of 

, M82. The two most compact radio objects in M82 have been investigated (41.95+575 

and 43.31+592). Using this recent epoch of data in comparison with our previous 
global VLBI observations and two earlier epochs of European VLBI Network observa- 

• tions we measure expansion velocities in the range of 1500-2000 km s _1 for 41.95+575, 
Qh and 9000-11000kms _1 for 43.31+592 using various independent methods. In each case 

the measured remnant expansion velocities are significantly larger than the canonical 
^ ■ expansion velocity (500 km s^ 1 ) of supernova remnants within M82 predicted from 

^ | theoretical models. 

In this paper we discuss the implications of these measured expansion velocities 
with respect to the high density environment that the SNR are expected to reside in 
within the centre of the M82 starburst. 

Key words: interstellar mediurn:supernova remnants - galaxies: individual :M82 - 

• galaxies:starburst - galaxiesunterstellar medium 



1 INTRODUCTION 

Studies of extra-galactic supernova remnants (SNR) are cur- 
rently limited by sensitivity, and to some extent angular 
resolution, nevertheless such studies provide unique insights 
into the early evolution of radio supernovae and supernova 
remnants. The extensiv e studies of remnants in our own 
galaxy (e.o. lGreerj200^1 are particularly valuable for testing 
details of the physical processes occurri ng in individual rem - 
nants (e. g. Braun, Gull & Perley 1987: lBrogan el aLl 2006). 
but they are limited by the fact that the youngest known 
galactic supernova remnant is over 300 years old and that 
the distances to galactic remnants, in many cases, remains 
quite unknown. 

Studies of extra-galactic SNR in starburst galaxies have 
a number of advantages. Firstly, the high star-formation rate 
results in a relatively large number of young supernova rem- 
nants with ages measured in decades, rather than centuries. 
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Secondly, as the starburst region is typically a kiloparsec in 
extent, and as distances to even the nearest starburst galax- 
ies are a few Mpc, the relative distances to each SNR will 
only vary by ~0.1 %. Finally, radio synthesis techniques en- 
sure that all the SNR, within an individual starburst galaxy, 
are observed with same angular resolution and sensitivity. 
Hence as all the SNR are essentially at the same distance 
this corresponds to a sample of SNR measured with constant 
linear resolution and surface luminosity limit. 

The star-formation rate (SFR) of a galaxy can be esti- 
mated using measurements in the ultra-violet . optical, infra- 
red and radio regimes (e.g. ICram et qZ.1ll998T> . However, as 
starburst galaxies usually contain large quantities of dust, 
the estimates of SFR from ultraviolet and optical data are 
uncertain due to dust extinction. M82 is no exception wit h 
A v typically 20-30 magnitudes (iMattila fc Meikle I l200ll) . 
Fortunately, radio methods are not affected by dust extinc- 
tion and hence the SFR can be estimated from the total 
1.4 GHz luminosity of the galaxy iCram et al. lll998l) . The 
total flux density of M82 at 1.4 GHz is ~8.5 Jy and using a 
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Table 1. Summary of VLBI observations of M82. 



Epoch Observation date Array Frequency Reference 



1 11 Dec 1986 (1986.95) Ef, Jb, Wb, Mc 1.4 GHz Pedlar et al, 1999 

2 02 Jun 1997 (1997.42) Ef, Jb, Mc, Nt, On, Wb, Cm, Tr 1.6 GHz Pedlar et al, 1999 

3 28 Nov 1998 (1998.91) VLBA, Yl, Ro, Go, Ef, Jb, Mc, Nt, On, Wb, Tr 1.6 GHz McDonald et al, 2001 

4 23 Feb 2001 (2001.15) VLBA, Yl, Ro, Ef, Jb, Mc, Nt, On, Wb, Tr 1.6 GHz This paper 
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Figure 1. Contour images of the compact radio source 41.95+575 from the global VLBI epochs observed on 28 Nov. 1998 and 23 Feb. 
2001. Both images have been convolved with a 3.3mas circular beam. The two epochs are contoured with linear multiples (—1, 1, 2, ... 
10) times 0.21mjybeam and 0.11 mjy beam - 1 for the epochs 3 and 4 respectively. 



distance of 3.2 Mpc, then the SFR(M>5M Q ) for M82 can 
be estimated using the relation given by Cram et al., to be 
~2.5 Mq yr _1 . Also the total FIR lumin osity implies a value 
of SFR(M> 5M<T)) close to 2M yr _1 dCram et al. I ll99St 
|Pedlarll200l . In addition the radio thermal free-free con- 
tribution is reasona bly well defined in M 82 and can be used 
to estimate a SFR JPedlar et al. 1119991) consistent with the 
above value. 

If the SFR is constant over ~ 10 7 years, then in princi- 
ple the type-II supernova rate can be derived directly from 
the SFR and the initial mass function (IMF), on the assump- 
tion that all stars with masses >8M become supernovae . 
Hence assuming a Miller-Scalo IMF iMiller fc Scalo. Ill979|) 
the supernova rate, i/ sn , is directly related to the SFR. This 
implies a v an ~0.08— 0.1 yr _1 in M82 and hence it is not sur- 
prising to find a relative ly large number of ra dio supernova 
remnants in this source ImuxIow et al. lfl994l) . 

Several methods have been used to estimate the v sn di- 
rectly from the radio parameters of the M82 supe rnova rem- 
nants iMuxlow et al. Ill994l : fPedlar et al. Ill999l) . The sim- 
plest method is to assume that all the M82 remnants which 
are brighter and smaller than the Cassiopeia A are younger 
than 330 years. Alternatively if a constant expansion veloc- 
ity is assumed (e.g. 5000-10000 km s " 1 ) this gives the age s 
of the remnants and hence the v „ „ iIMuxIow et al. lll994) . 
Finally Ivan Buren fc Greenhousel dl994l) estimated the v 3n 



directly from remnant luminosities. Each of these methods 
imply v an between 0.05 to O.lyr -1 , which is consistent with 
the SFRs discussed above. 

These supernova rates imply that the ^50 observed 
remnants have ages typically of a few hundred years. How- 
ever some doubts have been raised concerning the evolution 
of the remnants in M82. Kronberg, Biermann & Schwab 
(1985) claim, using statistical arguments, that in a large 
fraction of the remnants the radio luminosity is decaying 
at less than 0.1 % per year. If true this could suggest that 
the remnants are over a thousand years old . Also theoreti- 
cal studies bv lChevalier fc Fransson"! (1200 ll) have suggested 
that the high pressure of the interstellar medium in M82 
would result in the supernova remnants expanding at only 
500 km s -1 , hence as the observed sizes are typically a few 
parsecs, this again would imply ages of thousands of years. 

In order to investigate these problems it is essential that 
the expansion velocities and possible decelerations of the 
M82 remnants. This is possible using Very Long Baseline 
Interferometry (VLBI) techniques which enables their sub- 
parsec structure to be measured. Pedlar et al. (1999) used 
the European VLBI network (EVN) to produce the first 
VLBI images of the M82 SNR wit h 20 mas (0.3 pc) resolu- 
tion. This work was followed up bv lMcDonald et al. I (1200 ll) 
at even higher angular resolution (5 mas «0.1pc) using the 
Global VLBI network. 
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Figure 2. Contour images of the radio source 43.31+592 from the global VLBI epochs observed on 28 Nov. 1998 and 23 Feb. 2001. Both 
images have been convolved with the same 4 mas circular beam. The two images have been contoured with multiples of — 1, 1, 1.414, 2, 
2.828, 4 & 5.657 X O.lmjybeam - 1 and 0.08 mjy beam^ 1 respectively. 



In this paper we report on further global VLBI mea- 
surements taken over two years after the McDonald et al. 
measurements. They provide a 4th epoch of VLBI measure- 
ments and a 2nd epoch of global VLBI images of these SNR. 



2 OBSERVATIONS 

The latest global VLBI observations were made on the 23 
February 2001 using an array with 19 elements at a wave- 
length of 18 cm for a total of 18 hr. The array used con- 
sisted of the ten antennas of the Very Long Baseline Array 
(VLBA) in the USA, 7 antennas from the EVN, a single Very 
Large Array (VLA) dish and the NASA Deep Space Network 
(DSN) antenna at Robledo in Spain. This network was de- 
signed to be as close as possible to that u sed to make simi- 
lar ob servations on 28 November 1998 bv iMcDonald et al. I 
fcOOlT) which used a 20 element network. The NASA DSN 
antenna at Goldstone USA was not available for the observa- 
tions presented here, resulting in a reduction of the number 
of baselines to 171 compared to the 190 baselines used in 
the 1998 epoch. Due to the high declination of the target 
source, +69.7°, and long duration, 18 hr, of these observa- 
tions the non-availability of the DSN Goldstone antenna did 
not significantly compromise the excellent u-v coverage of 
this experiment. 

These data were taken in spectral line mode using a to- 
tal of 128 channels each with a bandwidth of 0.125 MHz. 
This yields a total bandwidth for the observations of 
16 MHz. 

Observations of several sources were made in addition 
to the main target source which was the central kiloparsec 
of M82. Throughout these observations 4.1 minute scans of 
M82 where interleaved with 2.1 minute scans of the nearby 
calibration source J0958+65 which was used for phase ref- 



erencing. During the observing run the bright calibration 
sources 3C 84 and J0927+39 were also observed and used 
in fringe fitting and bandpass calibration. The source 3C 84 
was also used in isolation as a flux density scale calibrator. 

These observations were correlated using the 
Mark IV/ VLBA system at the National Radio As- 
tronomy Observatory (NRAO) in Socorro, New Mexico. 
These data were then reduced using the Astronomical 
Image Processing System (aips) which is provided and 
maintained by the NR AO, following the s ame procedures 
as those detailed in lMcDonald et ail (120011) . 



3 RESULTS & ANALYSIS 

Although the location of all of the bright compact sources 
were imaged, as in our earlier studv lMcDonald et al. I (1200 ll) 
only two were observable with high signal-to-noise, due to 
the reduction in brightness sensitivity at high angular res- 
olution. Consequently the results presented here will focus 
upon the two most luminous and compact sources in M 82, 
41.95+575 and 43.31+592, and compare these results with 
the previous 1.6 a nd 1.4 GH z global VLBI and EVN observa- 
tions p resented in lMcDonald et al. I i200lT) and lPedlar et al~\ 
( 1999). The observational details of each of these epochs are 
summarised in Table 1. In the rest of this section we will 
discuss the structural evolution of the two detected compact 
radio sources 41.95+575 and 43.31+592. 

In Fig. 1 we present two images of 41.95+575, from 
bo th the 28 Nov. 1998 glo bal VLBI observations presented 
bv lMcDonald et al. I i200lT) and a matched resolution image 
created from this second epoch of global VLBI observations. 
Both of these i mages have b een created using a Briggs ro- 
bustness factor feriggslll995l) of —2 and convolved to an an- 
gular resolution of 3.3 mas. The measured rms noise of the 
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Figure 3. Contour images from all four VLBI epochs of the compact RSN 43.31+592. The date of observation of each epoch is shown 
at the top left-hand corner of each image. All four epochs have been convolved with a circular beamsize of 15 mas, to match the earlier 
EVN-only epochs in 1986 & 1997. All four images are contoured with -1, 1, 2, 3, 4, 5, 6, 7, 8, 9 and 10 times 0.35 mjy beam -1 . The 
grey-scale shown for three leftmost images is linear and ranges between 0.1 and 2.2mjybeam . The grey-scale shown on the 2001 
epoch is also linear but in the range of 0.1 and 1.7mjybeam — 1 . On the right-most images three crosses are plotted to show the relative 
Gaussian-fitted positions of the three components in epoch 1. 
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Figure 4. Two 7.5 mas wide slices along the major axis of the 
compact source 41.95+575. The bold black line shows the flux 
density versus position for the source in the 1998.91 global VLBI 
epoch whilst the lighter grey line is the equivalent slice across the 
2001.15 epoch. The dot-dash line represent the sum of the three- 
component Gaussian fit to this flux density, with the positions of 
the centroids of the two dominant Gaussian components, in each 
epoch, marked as extended ticks on the x-axis. 



new 2001 epoch is 36 /iJybeam -1 . The source 41.95+575 is 
the most compact of the sources in M82. This source was 
o nly m arginally resolv ed in the earlier EVN measurements 
of IPedlar et al. I ll 19991) , however both these new VLBI ob- 
servations and those of I McDonal d et al. I (|200l) clearly re- 
solve this source showing it to exhibit an elongated, bi-polar, 
structure unlike the shell-like structure commonly associated 
with SNR. The total flux density of 41.95+575 (including the 
diffuse emission surround the two compact components) has 
evolved between these two global VLBI epochs from 34.8+ 
0.1 mjy in Nov. 1998 to 28.9+0.3 mjy in Feb. 2001. 

The second most compact radio source within M82 is 
the shell-like remnant 43.31+592. Fig.|5|shows images from 
both of the global VLBI observations of 43.31+592. Unlike 



41.95+575, this source has a well defined shell structure and 
is much more typical of an SNe remnant. Although many of 
the features of the shell are common to both images, at 
the resolution of these global VLBI observations (few mil- 
liarcsec) the combination of structural evolution of the ring 
and image fidelity limitations resulting from the incomplete 
sampling of the u-v plane and non-linear effects of decon- 
volution, result in ambiguities when tracking individual fea- 
tures between the two epochs. Hence these images are also 
presented convolved with a 15 mas circular beam (Fig.|3J. 
At this lower angular resolution the effect of individual low 
signal-to-noise components is minimised; this also allows the 
two global VLBI observations to be directly compared with 
the early EVN-only observations made in 1986 and 1997 
JPedlar et al. Ill999n . 



3.1 The expansion of 41.95+575 

From our earlier measurements it was clear that the expan- 
sion velocity of this object was relatively low, and hence 
measuring the expansion velocity of this object over a two 
year baseline is particularly challenging. Furthermore the 
problem is exacerbated by changes in the small scale struc- 
ture which give rise to problems in finding common reference 
points between the two epochs. 

Several methods were used in an attempt to detect ex- 
pansion in 41.95+575, including taking slices along the ma- 
jor axis of the source as well as fitting the positions of the 
two peaks with 2-D Gaussian components. Figure^] shows a 
flux density distribution of a 7.5 mas wide, one-dimensional, 
slice along the major axis of 41.95+575. In this figure the 
lower, lighter grey line represents the latest global VLBI 
epoch (2001.15) and the upper solid dark black line the first 
global VLBI epoch (1998.91). Although the flux density evo- 
lution of this source between the two epochs is clearly shown 
in this figure, the small expansion of the two peaks is not 
easily discerned. However, a small increase in the separation 
is measured in the Gaussian fitted centroids of the two domi- 
nant components (the positions of these centroids are shown 
as extended ticks on the x-axis). This increased separation 
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Table 2. Fit of two Gaussians to the source 41.95+575 for each of the global VLBI epochs (epoch 3, 1998.91 and epoch 4, 2001.15). We 
have not included measurements derived from the earlier EVN-only observations in which this source is not resolved. The quoted flux 
densities, both peak and total are derived from the Gaussian fits of the two components and do not include the diffuse radio emission. 
The Right Ascensions and Declinations quoted are relative to 09 h 55 m , 69°40'(J2000). It should be noted that the absolute astrometric 
precision of these positions has been degraded due to the use of self-calibration, although the relative precision at each epoch has been 
preserved. 



Epoch 


Gaussian 
comp. 


Size 
(mas) 


PA 

(deg) 


RA 
(J2000) 


Dec 
(J2000) 


Speak 

(mjy/bm) 


Stot 

(mjy) 


Distance Between Peaks 
(mas) (pc) 


3 


1 
2 


12.7x7.6 
16.2x6.6 


23±2 
21+1 


50.68813±0.00002 
50.68529±0.00002 


43.7764±0.0001 
43.7682±0.0002 


1.80±0.05 
1.38±0.04 


17.9+0.5 
15.2+0.5 


16.6+0.2 


0.22+0.03 


4 


1 
2 


15.7x11.3 
20.2x11.2 


47+3 
36±2 


50.68789±0.00002 
50.68502±0.00004 


43.7672±0.0001 
43.7583±0.0002 


1.14±0.02 
0.71±0.02 


12.9±0.4 
9.1+0.5 


17.7+0.2 


0.24+0.03 



Table 3. Derived radii of SNR shell 43.31+592. The first two columns list radii of peak and 50% level of the flux density when averaged 
within concentric annuli. The flux density versus radius distribution for each epoch is graphically shown in Fig 5. The final column shows 
radii measured for each epoch derived using the analytical 'Common Point Method' which is described in more detail in Section 3.3.1. 
In each of these cases the values listed have been derived from the matched, 15 mas angular resolution images of 43.31+592. 





Annuli (at peak) 


Annuli (at 50% of peak) 




CPM 


Epoch 


Derived radius 


Derived radius 


Derived radius 


Derived radius 


Derived radius Derived radius 




(mas) 


(PC) 


(mas) 


(PC) 


(mas) 


(PC) 


1 


13.1 + 1.0 


0.179 ± 0.013 


22.2 ± 1.0 


0.296 ± 0.013 


21.0 + 0.4 


0.326 ± 0.006 


2 


19.9 ± 1.0 


0.265 ± 0.013 


29.8 ± 1.0 


0.397 + 0.013 


29.4 + 0.6 


0.456 ± 0.009 


3 


20.1 ± 1.0 


0.268 + 0.013 


30.5 ± 1.0 


0.407 + 0.013 


31.9 + 0.6 


0.495 ± 0.009 


4 


20.4 + 1.0 


0.272 ± 0.013 


32.4 ± 1.0 


0.432 ± 0.013 


33.1 + 0.6 


0.514 + 0.009 



Gaussian fits to the two brightest components within the 
image plane. The results of this fitting are shown in Ta- 
ble|21 The expansion derived from these fitted positions is 
1.1+0.3 mas over the 2.24 year baseline implying a radial 
expansion velocity of 1500+400 kms -1 . 

The expa nsion along the major axis has been meas ured 
in the past bv lTrotmanl Jl99fJ).lMcDonald et al.\ J200ll) and 
more recently in iMuxlow et al. I fe005bl) ~ The results pre- 
sented here are consistent with these results. It is necessary 
to note that as this source appears to be bipolar in nature 
the expansion along the major axis may be greater than the 
value quoted here due to orientation effects. Consequently, 
and considering the large errors in measuring the expansion 
of 41.95+575 over only a 2.24 year baseline, we tentatively 
conclude that this source is radially expanding at a rate of 
~ 1500-2000 kms" 1 . 



3.2 The expansion of 43.31+592 

A comparis on between the 198 6 and 1997 EVN images of 
43.31+592 (IPedlar et al. 1119991) showed an obvious increase 
in size consistent with an expansion velocity of ~ 10000 
kms -1 . The continuing expansion of this source is evident 
in Fig. 3. 

More detailed expansion measurements of this source 
using the full resolution Global VLBI data poses problems 
due to evolution of the small-scale structure which did not 
allow accur ate Gaussian fitting to th e four knots of emission 
identified in iMcDonald et al. M200ll) . Instead, the centre of 
the source in the maps was found for each epoch and the flux 
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Figure 5. Radial flux density profile of the SNR 43.31+592 for 
each epoch of observation. The left-hand panel shows four pro- 
files, for each epoch, derived from images convolved with a 15 mas 
synthesised beam, whilst the right-hand panel shows similar pro- 
files for the two global VLBI images, 1998.95 and 2001.15, when 
convolved with a 4 mas circular beam. 



is measured to be 0.7+0.2 mas year - , which is equivalent to 
a projected radial expansion rate of 2050+500 kms -1 . 

In addition to fitting the one-dimensional flux den- 
sity distribution versus distance along the major axis of 
41.95+575, we have performed extensive two-dimensional 



6 Beswick et al. 

was integrated in a series of 1 mas-thick annuli at increasing 
radii from the centre. The geometric centre of 43.31+592 was 
established in each image as the position which was equidis- 
tant from Gaussian fitted components within the ring. A 
plot of the integrated flux in each annulus against the dis- 
tance from the centre is plotted in Fig.|^] for all epochs at 
a matched angular resolution of 15 mas and for the two lat- 
est epochs of global VLBI observations at a resolution of 
4 mas. The radius of the flux density maximum and 50% of 
this value for each of the 4 epochs, measured at a matched 
angular resolution of 15 mas, are listed in Tabled The ro- 
bustness of these derived radii was tested with respect to the 
chosen centre of the remnant by repeating this exercise for 
several centre positions up to a few beam-widths from the 
calculated centre. It was found that, although the radius 
of the peak of the radial flux density profile significantly 
changed and smeared when the nominated centre was more 
than a beam-width from the true image centre, the radius as 
measured at a point between ~30-70% of the peak flux did 
not vary significantly provided the central position deviate 
by more than a beam-width from the true remnant centre. 

As can be seen in both Tableland the left-hand panel 
of Fig.|21 a distinct expansion is observed in the peak ra- 
dial flux density distribution between the 1986.95 epoch and 
latter epochs in the lower (15 mas) resolution data. How- 
ever little discernible expansion of the peak flux density is 
recorded over the shorter time lines between the three lat- 
ter epochs, which is not unexpected since an assumed ex- 
pansion velocity of 10000 kms" 1 would imply an increase 
in the peak radius of ~2.5mas over the 3.7 yrs between 
1997.42 and 2001.15. Using the higher 4 mas angular resolu- 
tion global VLBI data, the measured peak of the radial flux 
density distribution (Fig0 in 1998.91 is 19.5 mas compared 
to 22.5 mas in 2001.15. This implies an angular expansion 
speed of 1.34+0.32 mas yr" 1 , corresponding to a radial ex- 
pansion speed of 7800±1900 kms" 1 . 

An alternative measure of the radial flux density dis- 
tribution is the radius at a percentage of the peak flux. In 
Tabled in addition to the peak flux density radius, we quote 
the half peak flux density radius for each epoch of observa- 
tions as measured from the lower 15 mas resolution images. 
Whilst neither of these values is a true representation of 
the radius of the SNR's radio emission, both are directly 
proportional to this value and can be used to derive the 
expansion velocity of the SNR (this is discussed in further 
detail in the following section). In Fig.|S|radial expansion ve- 
locities between each epoch of observations have been plot- 
ted, along with an independent measurement of the velocity 
of 43.31+592 made by comparing deep Multi-Element Ra- 
dio Linked Interferometer Network (MERLIN) 5 GHz ob - 
servations made in 1992 and 2002 (|Muxlow et al. Il2005ah . 
The mean radial expansion velocity of these VLBI measure- 
ments is 9025+380 kms" 1 compared to 8750+400 kms" 1 , 
recorded over a similar period using MERLI N. These result s 
are consistent with the results i n |M^Donak4rf al. I i200ll) , 
IPedlar et al\ fjl999l) and lMuxlow et al. I fj2005aTl . 

3.3 The expansion of 43.31+592 measured using 
an analytical method 

As well as using the conventional methods to measure shell 
sizes, we have also applied a novel analytical technique de- 
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Figure 6. Derived velocities for the SNR 43.31+592. Horizon- 
tal bar represents the time baseline between size measurement, 
vertical bar is the error in measuring the shell radius at 50% 
peak flux density value. Velocities measured between all available 
epochs have been included. An additional independent measure 
(g) of the velocity of 43.31+592 by Muxlow et al, (2005), using 
MERLIN, is also included for completeness. 



veloped to parameter ise the expansion of SN1993J in M81 
IMarcaide et al\ fcOOfJ) . 

When trying to analyse the expansion of a shell-like 
object it is important to take into account the continuously 
increasing ratio between the size of such a shell and the (con- 
stant) beam used for the convolution of the CLEAN compo- 
nents to obtain the final image. As the source is expanding, 
the ratio of the convolved beam size and the source size will 
become proportionally smaller, this fact can produce biases 
in the determination of the expansion parameters. Whilst 
methods for the correction of this beam-re lated bias in the 
imagin g analysis have been developed by IMarcaide et al. I 
(1997), such solutions tend to produce other undesirable 
new biases in the data. However, a new method for the 
determination of the expansion of radio-emitting shell-like 
objects that so l ves th is beam bias has been developed by 
IMarcaide et al\ J2006I) . When applied to shell-like objects, 
the method has been shown to work in a way totally inde- 
pendent of the a priori beam used. This method, called the 
'Common Point Method' (CPM), is based on mathematical 
properties of the convolution of angular symmetric images 
with Gaussian components. T his method is described in de- 
tail in IMarcaide et al\ ||2006). 

In essence, the CPM exploits the fact that there is a 
point in the angular average of such circularly symmetric 
images that is stable with respect to changes in the con- 
volving beam (see Fig.|7J. This stability can be used in an 
iterative process in order to obtain a magnitude that will 
be perfectly proportional to the radius of the source, with 
a proportionality constant close to unity and will remain 
constant provided the expansion is nearly self-similar. 

The results from applying the CPM to these 43.31+592 
data are listed in Tableland are in good agreement with 
the other expansion measurements presented in this pa- 
per. The measured radial expansion velocity of the remnant 
43.31+592 using the CPM is 11000 ± 1000 kms" 1 , fitting to 
each of the four epochs and assuming a free expansion model 
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Figure 7. Qualitative representation of the properties on which 
the CPM is based. The solid line represents the radial profile of 
a given shell-like emitter, whilst the dashed lines show profiles 
of the angular average of the convolution of this shell with two 
different beams, one of which is 50% wider than the other. The 
outer point where the dashed curves coincide is very stable under 
beam changes and, thus, can be found to be related with the 
radius of the source within a given iterative process. Note also 
that the maxima of both dashed curves is at a radius smaller 
than the true source radius. 




Figure 8. Measured radius of 43.31+592, using the CPM, for 
the four epochs. We also show five solid lines representing fitted 
models with fixed expansion parameters ranging from 0.4 to 1.0 
(free expansion). The vertical dashed line marks the epoch of the 
first detection of this source in 1972 by Kronberg & Wilkinson 
(1975). 



{i.e. linear expansion). Hence if the expansion is linear, the 
supernova explosion took place in ~1962. 

It is possible to fit a wide range of deceleration parame- 
ters (5) to these four derived shell r adii, since the birth date 
of t his source is unknown . Following lMcDonald et aC\ (|200ll ) 
and lHuang et al. I al994T) we parameterise the size evolution 
of the remnant as, D — kT 5 , where D and T are the size 
and age of the remnant respectively, k is a constant and S is 
the deceleration parameter. In Fig. [H] the derived shell radii 
are plotted with various fitted parameterised deceleration 
curves overlaid. It should be noted that in 1972 (marked by 
the vertical dashed line in Fig.El, 43.31+592 was dete cted in 
radio images made bv iKronberg fc Wilkinson! lll975l) . This 
early radio detection of this source implies a lower limit on 
S of 0.68. 



4 DISCUSSION 

From these global VLBI measurements a radial expansion 
velocity of 1500+400 kms™ 1 is derived for the most compact 
radio source in M82, 41.95+576. This source is, however, 
rather anomalous on account of both its high radio luminos- 
ity (compared to the other supernova remnants within M82) , 
rapid flux density decay and highly elongated, bipolar, radio 
structure. In particular the radio structure of 41.95+576 is 
highly atypical for a SNR. As yet it is uncertain what this 
source is, however it is clear that it is inappropriate to treat 
it as a typical SNR. Th e nature of 41.95+576 is discussed in 
iMuxlow et oTI <2005bl) . 

The well-defined shell structure of 43.31+592 is more 
typical of SNR and radio supernovae. As can be seen in 
Fig.|]3|th.e present observations are consistent with near-free 
expansion of the shell remnant 43.31+592 with a veloc- 
ity close to 10000 km s _ . This expansion velocity is more 
than an order of magn i tude larger than that predicted by 
IChevalier fc Fransson I l|200lft which was largely a conse- 
quence of the high density and pressure assumed for the M82 



ISM. The minimum deceleration parameter that is consis- 
tent with this expansi on and the first known record ing of the 
radio source in 1972 iKronberg fc Wilkinso"r]ll975l) is 0.68, 
which is close to the J50.73 limit placed on this paramete r 
from earlier VLBI observations bv lMcDonald et al. I i200ll) . 

4.1 Interstellar medium pressures in M82 

IChevalier fc Fransson I ^2001) assume that the pressure in 
the starburst region of M82 is 10 7 cm -3 K and from this 
predict low SNR expansion velocities of ~500kms~ . While 
there is no doubt that part of the starburst must b e at th is 
pressure (e. g. compact Hll regions iMcDonald et al. H200ll) ). 
it seems likely that much of the starburst is not in static pres- 
sure equilibrium and that a range of pressures are present. 

As the photoionised component in M82 has a tem- 
perature of ~ 10 4 K, the pressure assumed by Chevalier 
& Fransson would require that the density of this gas be 
~10 3 cm -3 . However, observations of free-free absorption 
against individual r emnants in M82 gi ve emission measures 
of ~5x 10 5 pc cm -6 JWills et al. 111997!) . Consequently if the 
gas number densities were to be ~10 3 cm -3 the observed 
emission measure would imply ionised gas path lengths of 
only ^0.5 pc. Given that the extent of the starburst region 
is at least 500 pc, unless the ionised gas has contrived to 
occupy only a thousandth of the volume of the starburst in 
the immediate vicinity of each remnant, it is difficult to see 
how the ionised component can provide sufficient pressure 
to slow the SNR expansion to 500 km s -1 . 

Chevalier & Fransson suggest that the SNR may be 
confined by the interclump medium of molecular clouds in 
M82 and note that many of the remnants are in the line 
of sight to regions of strong CO emission. This is, however, 
not conclusive evidence that the remnants are embedded 
within the clouds, and often the Hi absorption spectra mea- 
sured directly against the remnants JWillsll998tfWills et al. I 
1998) do not support this interpretation (see discussion of 
43.31+592 below). Furthermore, studies by Weifi, Walter 
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& Neininger (1999) have shown that the molecular gas in 
the starforming regions of M82 has a kinetic temperature of 
— 150K and number densities of ~ 10 3 cm -3 which corre- 
sponds to a pressure two orders of magnitude less than that 
assumed by Chevalier & Fransson. 

Clearly the ISM of M82 is complex, and to assign a 
single pressure to the starburst region may be unrealistic. 
Hence it is likely that most of the SNR in M82 are em- 
bedded in regions with pressures significantly lower than 
10 7 cm -3 K, especially since pre-supernovae stellar winds 
from the massive progenitor stars of these remnants are 
likely to have driven away much of the material within the 
vicinity of these young supernova remnants prior to the su- 
pernova phase beginning. Consequently, it is highly plausi- 
ble that the observed remnants are expanding into regions 
of the ISM with lower density than the average density and 
hence have velocities significantly in excess of 500 kms" 1 . 



4.2 Is 43.31+592 in low density region? 

The size, radio luminosity and lack of rapid flux variability of 
the source 43.31+592 appears to be typical of young SNRs. 
In addition it shows a well-defined shell structure (Fig.|2HH]l. 
As the shell seems to be in, or close to, free expansion (see 
Section 3.3) a simple constraint of < 2000 atoms cm -3 can 
be set to the external density by assuming that the mass 
of gas swept up is less than the mass ejected. The expan- 
sion velocity of this source clearly exceeds the 500 kms" 1 as 
discussed above. 

This possibility is sup ported by the lack of a low fre- 
quency turnover at 73 cm JWills et al. Ill997l) in the spec- 
trum of 43.31+592 even though it is in the line of sight 
to regions with strong free-free emission. This is consistent 
with the remnant being in front of, and not embedded in, 
the relatively dense ionised gas in the central region of M82. 
Additionally, even though relatively strong Hi abso rption 
(5 x 10 21 atoms cm" 2 ) is seen against this remnant dWillsl 
1998), the absorption has a relatively narrow width and the 
systemic velocity is consistent with absorption by a disk of 
Hi external to the central region of M82. Hence this could be 
further evidence that this source is not embedded in neutral 
gas in the central region of the starburst. Thus if this SNR 
is in a particularly low pressure region, it might possibly 
reconcile the difference between the theoretically predicted 
and measured expansion velocities. 

However it seems unlikely that such a scenario can be 
contrived to account for high expansion velocities in more 
than a few of the M82 SNR. Recent MERL IN observations 
jMuxlow et al. Il2005al: iFenech et al. 1120051) have indepen- 
dently confirmed the expansion velocity for 43.31+592 in 
addition to expansion velocities of several more remnants. 
Of these other remnants their observed radial expansion ve- 
locities range from -2000 kms" 1 to over 10000 kms , im- 
plying that neither the relatively slow expansion velocity of 
41.95+575 nor the much higher velocity of 43.31+592 are 
particularly atypical. 

Hence, we see no reason to revise the supernova rates 
or remnant ages that we have inferred previously. However, 
we can offer no simple explan ation for the lack of var iability 
seen in the compact sources jKronberg et al. Ill985ft . Over 
the next decade any variability will be constrained by new 



observations using the next generation of high sensitivity 
radio interferometers, such as the EVLA and e-MERLIN. 



5 CONCLUSIONS 

Using these global VLBI observations, the two most compact 
radio components within the centre of the M82 starburst 
have been imaged with angular resolutions of a few milliarc- 
seconds. This latest epoch of observations, in conjunction 
with our previous VLBI observations of M82, complete a 
~15 year timeline of VLBI observations of these compact 
sources. 

The most compact radio source in M82, 41.95+575, has 
been imaged using both these latest data and those pre- 
viously published in Mc Donald et al."] {2001). This compact 
radio source has an elongated, bipolar structure that is atyp- 
ical of a SNR. Between these two global VLBI observations 
separated by 2.2 year, this source has grown with a radial 
expansion velocity along its major axis and in the plane of 
the sky of 1500+400 kms" 1 . Between these two global VLBI 
observations the total flux density (including both of the two 
compact components and the surrounding diffuse emission) 
of 41.95+575 has reduced at a rate of 7.1+0.5% per year, 
slight ly lower t han the previously reported 8.5% per year 
jTrotmanlll99rjl . 

The shell-like remnant 43.31+592 has been imaged 
with synthesised beams of 4 and 15 mas, affording com- 
pa risons with both t he pre vious global VLBI observations 
of iMcDonald et al. I feOOll) and t he earlier EVN-only ob- 
servations of |Pe*dlareFaIn l|l999fi . Various methods have 
been used to derive the expansion of this radio shell, all of 
which provide consistent radial velocities of ~ 10000 kms" 1 
between each of the four independent observations. This 
value is consistent with expansion measurements of this 
source from observatio ns using other radio interferometers 
jMuxlow et al. Il2005al) . 

Over the ~15 year timeline of these observations 
43.31+592 appears to be evolving at a rate close to 
that expected for free-expansion. Using the earliest def- 
inite radio detection of th is compact source in 1972 
jKronberg fc Wilkinson! Il975l) these observations imply a 
lower limit on the deceleration parameter of 0.68. Contin- 
ued VLBI observations over an increased time baseline will 
determine the deceleration rate of this source, as well as al- 
lowing the interaction of the SNR shell with its surrounding 
ISM to be further investigated. 

Higher sensitivity observations using e-MERLIN and 
broadband VLBI, combined with an increasing time base- 
line, will enable the expansion velocities of many more rem- 
nants in M82 to be determined and hence constrain both the 
supernova remnant parameters and the properties of the in- 
terstellar medium in M82. 
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